Double parton scattering as a source of quarkonia pairs in LHCb 
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Recent results on production of J/ip pairs in LHCb initiated discussion on double parton scattering 
(DPS) contribution to this process, which is not yet elaborated well. In this short report J/ipXc an( A 
J/ipT modes are proposed for DPS studies. Estimations for TT and J/tpT production in LHCb are 
presented. 
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Recent observation of J/ip-meson pairs at LHCb stimulated discussion on phenomena contributing this process. 
On the one hand, there is leading order (LO) calculation of the gg — > J/ipJ/ip process in the color-singlet (CS) model 
[3-[3]- It leads to the total cross section at 7 TeV energy of 24 nb and to 4 nb in LHCb conditions. On the other hand, 
high flux of incoming partons at LHC energies leads to significant probability of scattering of more than one parton 
' pairs in the same proton-proton collision [5]. Estimations of this double parton scattering (DPS) contribution to the 
J/ip pairs production lead to the cross section of approximately 2 nb in the LHCb acceptance 0]. Both SPS and 
DPS estimations are of the order of first experimental measurement in LHCb Q, a(pp — > J/ip J/ip + X) = 5.6 ± 1.1 
nb. Meanwhile, both these predictions have some uncertainties. In SPS calculations they are induced mainly by 
dependencies on a s and m c parameters. DPS estimations does not account for correlations in double parton density 
functions and depend on a phenomenological parameter er cff . Thus further investigation is desirable. 

In this work other di-quarkonia final states are considered in addition to the J/ip J/ip one. It will be shown that for 
J/ipXc and J/ipY production DPS should be main source at least at low-invariant-mass region. 

Assuming factorization of two hard partonic processes A and B one can write inclusive cross section of a double 
parton scattering process in a hadron collision in the following form Q : 
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where r^(xi, x 2 ; h»i, h 2 ; Q 2 , Q 2 ) is double parton distribution function (PDF), crf^ B — partonic cross sections of 
processes in question, b — impact parameter and m = 2 for different partinic subprocesses and 1 — for identical. It 



,. is usually assumed that longitudinal and transverse components of the PDF can be decomposed in the following way: 

T^x^hiMiQlQl) = ^(x 1 ,x 2 ;Qig2)/(b 1 )/(b 2 ), (2) 
and longitudinal component D % ? (x\, x 2 ; Q 2 , Q 2 ) is taken as a product of two independent single parton distributions, 
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D^(x u x 2 ;QlQ 2 ) = D l h ( Xl ;Q 2 )Di(x 2 ;Q 2 ). (3) 
This leads to a well-known simple expression for DPS cross-section in which no di-parton correlations are involved: 
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CDF and DO measurements give a c s = 14.5 mb, which is roughly 30% of the non-single diffraction (NSD) cross 
section at the Tevatron (» 48mb). The NSD cross section at LHC is only slightly higher (« 51mb). This supports an 
assumption that <j e g weakly depends on the total energy of interaction [l(| • Nonetheless dependence on the resolution 
scale and consequently on the partonic process can be significant. Model of multiple interactions implemented in 
Pythia8 [ll[ MC generator tries to accounts for this dependence and predicts a c g for double charmonia production 
of 30 mb (cr c ff = ctnsd/ (/impact))- Further we will use usual value of 14.5 mb. 

LO calculations in the CS model lead do the value of double prompt J/ip production cross section in LHCb 
acceptance (2 < y < 4.5) of approximately 
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4 nb. (5) 
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FIG. 1: Feynman diagrams for gg — >• qqqq process. 



This result depends on the ipj/^^O), rn c and a s parameters. First of them, ^jm(O), can be determined rather precisely 
from the leptonic width. m c was set equal to mj/^,/2. LO running a s at the transverse mass scale was used. Taking 
experimentally measured cross section of single J/tp production in LHCb of 10^b [l2| and using expression (j4| one 
gets for the DPS contribution approximately 

^j/iPJ/tP+x = 2nh (6) 

Sum of these contributions agrees well with the experimentally measured value js] 

a PP^ J /^ J ^+ x = 5.6±l.lnb. (7) 

Several methods to distinguish DPS contribution from the SPS one were proposed [H, 0] . However it would be easer 
to deal with process in which one of these contributions in suppressed compared to another. 

Crucial feature of the CS model is presence of so-called selection rules. According to the C-parity conservation, 
C-parity of the final state must be the same as those of 2 initial gluons, which is C-even as they are in CS combination. 
That is why production of J/ipXc in SPS is expected to be suppressed. DPS cross section has no suppression in this 
mode and should be significant as it is known that about 50% of single J / ip mesons originate from the feeddown from 
the Xc decays. Thus observation of J/V'-pairs accompanied by a photon (from the Xc J/^Pl decay) would be a 
signal of DPS contribution. 

Let us consider T(lS')-meson pair production. Calculations analogous to J/^-pair production in [2j lead to the 
total cross section at 7 TeV of 31 pb and in LHCb acceptance: 

a^ rr+x = 8.7 pb, (8) 

while estimation of DPS cross section based on the experimental data on single T production at LHC [l3[ gives 

a«& TT +* = 0.4 pb. (9) 

One sees that in the T(15)T(15) mode SPS should dominate. Also feeddown from T(2S) and T(3S I ) decays should 
increase SPS prediction while it is already accounted for in the DPS estimation. 

There are basically 3 types of Feynman diagrams for the gg — > qqqq process (see fig. [1]). Diagrams of first type 
contain 1 fcrmion loop while diagrams of second type contain 2 fermion loops connected by an intermediate gluon. 
In the diagrams of second type final particles are coupled with 2 gluons and consequently these diagrams do not 
contribute to 5-wave pairs production due to the C-parity conservation. Diagrams of the third type correspond to 
gluon fragmentation and contribute to CO-states production only. Thereby only diagrams of first type contribute 
CS J /if) or T pair produced. But there are no LO diagrams contributing CS J/ipY combined production. This final 
state is however accessible through XcXb production followed by Xc — > J I ip + X and Xb — >• Y + X decays. Meanwhile 
XcXb production does not exhibit kinematical peak near the threshold and is not expected to be significant in the low 
invariant mass region. Prediction of J/i/jY production in CO mechanism made in [4| give in LHCb acceptance 

^j/tf.T.+x = 2 pb; (1Q) 
while DPS contribution found using expression (jj) gives approximately 
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TJ/^+X = 75pk (u) 



So in the J/ipY mode prediction for DPS exceeds significantly those for SPS one. It is interesting to notice that for 
SPS crg^ T < cr|p T s < crg^ while for DPS cr£J s < < . 

Thereby cross section of J/ipJ/ip production is saturated by both SPS and DPS contributions, while for TT and 
J I tpT final states only one of these regimes prevail. 
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